Dynamical polarizability of graphene beyond the Dirac cone approximation 
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We compute the dynamical polarizability of graphene beyond the usual Dirac cone approximation, 
integrating over the full Brillouin zone. We find deviations at hoj = 2t (t the hopping parameter) 
which amount to a logarithmic singularity due to the van Hove singularity and derive an approximate 
analytical expression. Also at low energies, we find deviations from the results obtained from the 
Dirac cone approximation which manifest themselves in a peak spitting at arbitrary direction of the 
incoming wave vector q. Consequences for the plasmon spectrum are discussed. 

PACS numbers: 63.20.-e, 73.20.Mf, 73.21.-b 



I. INTRODUCTION 

Graphene is a novel two-dimensional system with many 
outstanding mechanical and electronic properties 1 . Espe- 
cially the early observation of the ambipolar held effect 2 
and of the odd integer quantum Hall effect 3 have stim- 
ulated enormous research on the electronic structure of 
graphene. Only recently, the fractional quantum Hall ef- 
fect was seen in suspended graphene 4 . For a review of 
this newly emerging branch of condensed matter physics, 
see Ref. 5 . 

To understand the unusual electronic properties of 
graphene, it often suffices to discuss the charge suscep- 
tibility. The static polarizability at hp, e.g., gives the 
Thomas-Fermi screening length, important for transport 
properties 6-8 whereas the dynamical polarizability at 
zero-wave number can explain the phonon softening 9 at 
the T-point. It is also used for the understanding of struc- 
tural inhomogeneities in graphene, so-called ripples 10 and 
the van der Waals interaction between graphene layers 11 . 

For neutral graphene, the polarizability at zero tem- 
perature was first calculated by Gonzalez et al. 12 , the ef- 
fect of temperature was discussed by Vafek 13 and vertex 
corrections were considered in Ref. 14 . For a gated system 
with finite chemical potential, the first expressions were 
given by Shung 15 in the context of graphite and later 
by Wunsch et al. 16 and Hwang and Das Sarma 17 . Also 
the extension to finite temperature has been performed, 
even though a closed analytical expression is then - as in 
the neutral case - not possible, anymore 18 . Recently, the 
polarizability was discuss in the presence of a magnetic 
field 19 and gapped graphene 20,21 . 

All these results originate from the Dirac cone approx- 
imation in which the energy dispersion of the hexagonal 
lattice is linearized around one of the two Dirac points 
where the valence and conduction band touch. But cor- 
rections to this approximation have to be included to dis- 
cuss e.g. the recently measured absorption of suspended 
graphene in the visible-optics regime 22 , which is related 
to polarizability via the continuity equation. This has 
been done in a perturbative treatment in Ref. 23 . The 
optical properties of graphite were calculated in Ref. ? . 



Here, we want to extend the previous calculations to 
the full Brillouin zone of the hexagonal tight-binding 
model. We certainly expect deviations at large energies 
where the Dirac-approximation does not hold anymore. 
But the main purpose is to test whether the diverging 
density of states at the M-point (van Hove singularity) 
leads to consequences on the collective excitations of this 
system. 

Our interest is motivated by the recent findings of 
an additional plasmon mode that emerges at around 
4.7eV with a linear dispersion 25 which was observed 
on freestanding graphene by electron energy loss 
spectroscopy 26 . A first guess is to associate this mode 
to the van Hove singularity which in the charge suscep- 
tibility shows up at 2t 5.4eV, t denoting the tight- 
binding hopping parameter. Including excitonic effects, 
the prominent absorption peak shifts to 4.5eV 27 , thus 
suggesting that the van Hove singularity is indeed the 
origin for this new plasmon mode. 

Apart from the high-energy corrections stemming from 
interband transitions, we also look at the intraband con- 
tribution at low energies and find that even there devia- 
tions from the Dirac cone approximation occur. By this, 
we complement a recent work, where the plasmon dis- 
persion is discussed by also considering the full Brillouin 
zone 28 . 

The paper is organized as follows. In Sec. II, we intro- 
duce the model and notation and define the polarizability 
of graphene. In Sec. Ill, we discuss the imaginary part 
of the polarizability which will only involve one numeri- 
cal integration. We first treat the interband contribution 
where we especially focus on the behavior around the M- 
point where the van Hove singularity occurs. We then 
discuss the intraband contribution and the different be- 
havior at certain directions of the incoming wave vector 
q. In Sec. IV, we obtain the real part of the polarizability 
via the Kramers-Kronig relation and discuss implications 
on the modified plasmon spectrum due to the inclusion 
of the whole Brillouin zone. We close with conclusions 
in Sec. V and give details on the analytical evaluation of 
the polarizability around the M-point in an appendix. 
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II. THE EFFECTIVE MODEL AND THE 
POL ARIZ ABILITY OF GRAPHENE 

The Hamiltonian of a hexagonal graphene sheet in 
Bloch spinor representation is given by 



a) 



H = £(-tff k - ^2x2) , Hk = 4 t) k ) ^ ' 



(1) 



with t w 2.7eV the tight-binding hopping parameter, 
the chemical potential and Vk = (ak,frk) T , «k and 6k 
being the destruction operators of the Bloch states of 
the two triangular sublattices, respectively. Further, we 
have 0k = ^2s e~^ <5i ~ <53 )' k , where Si denote the three 
nearest neighbor vectors. Here we choose them to be 

5i = |(-l,\/3), 5 2 = |(-1,-V3) , «5 3 = a(l,0), (2) 

with a ~ 1.42Abeing the nearest carbon-carbon distance. 
The Brillouin zone is then defined by the two vectors 
bi = 27r/(3a)(l,\/3) and b 2 = 27r/(3a)(l, — y/3), see Fig. 
1 a). The dimensionless eigenenergies thus read 

|</> k | = \Ji + 2cos(v / 3fc 2/ a) + 4cos(V3k y a/2) cos(3fca;a/2) . 

(3) 

In terms of the bosonic Matsubara frequencies ui n — 
2nn/(3 {(3 = l/fc^T, ft = 1), the polarizability in first 
order is defined as 



pW(q,^„) = i jf dre^ T (p(q,r)p(-q,0)) 



(4) 



where A denotes the area of the graphene sample and 
the density operator is defined as pq = /Oq + e 4q ' <53 ( o q with 

Pq = Ek,a c k, ff c k+q,<T (c = a, b). 

Hence, we obtain the general expression for the polar- 
izability 



P«(q,iu; n ) = — ^ f d 2 k V / M ,(k,q) 



(5) 



S.S 7 — ± 

n F (P*(k)) -n F (P s '(k + q)) 
P s (k) -P s '(k + q) +iw„ ' 

with E ± (k) = ±t\(j>k \ — fj, and n F (E) the Fermi function. 
For the band-overlap, we have 



/±(k,q) = -(l±Re 



s iq <53 . 



V k+q 



I0k| |0k+q| 



(6) 



Note that since we are summing over the entire Brillouin 
zone, only the spin-degeneracy g s — 2 has to be taken 
into account. 

For neutral graphene, fi = 0, there is no intraband 
contribution due to the canceling Fermi functions in the 
numerator of Eq. (5). Due to /_(k, q — > 0) — > 0, we 
further expect no interband contribution for q = 0. We 
finally note that for high energies uj > t, the phase fac- 
tor between the particle densities of the two sublattices, 
e »q-<5 3; i s crucial even in the long- wavelength limit q — > 0. 
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FIG. 1: a) The hexagonal and rhombical Brillouin zone b) The 
symmetrized rhombical Brillouin zone and its segmentation. 
The inner square refers to j = — , the outer triangles refer to 
j = +; additionally the values of s and s' are given. 



III. IMAGINARY PART OF THE 
POLARIZABILITY 

With the substitution ius n — > uj + iO, the imaginary 
part of the retarded susceptibility is written in terms of 
a delta function in the usual way. Determining the zeros 
of the argument of the delta function allows to perform 
the integration over k x analytically. The subsequent in- 
tegration over k y is then done numerically. We have also 
performed the direct summation of Eq. (5) of a finite 
system to check our results. 

The integration over the Brillouin zone can be split up 
into separate parts with slight modifications of the inte- 
grand (see below and Fig. lb)). The final domain is then 
given by < 3k x a/2 < tt/2 and < VSk y a/2 < ir/2 and 
the substitution x — sin(3fc x a/2) and y = sin(^/3k y a/2) 
can be performed. The resulting expression (see Eq. (8)) 
explicitly displays the inversion symmetry of q with re- 
spect to the q x - and q y -axis. The polarizability pW (q, u>) 
is also invariant under rotation of 7r/3, displaying the un- 
derlying lattice symmetry. We thus find the following 
symmetry: 

P«(|q|,7r/6 + <p,w) = P«(|q|,7r/6 - <p,u) (7) 

The subsequent plots thus only show four representative 
curves with < <f < tt/6. 



A. Interband transitions 

We shall first discuss the contribution of the interband 
transitions to the imaginary part of the polarizability. 
As explained above, we first perform the integral over 
k x — > x, thus eliminating the delta function. For neutral 
graphene, fj, = 0, at zero temperature T = this yields 
the following expression: 



(1) 2sgnHW3 / 2 \ 2 f 1 dy 



v 1 -y 

x Fj{ x iiU'i s Qx, s'q v ) 

j=± s,s'=± 0< Xl <l \A ~~ x i Vi s Qx,s'qy)\xi 



(8) 
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FIG. 2: (color online): Left hand side: The imaginary part of 
the polarizability ImPW(|q|, tp, w) as function of the energy 
ui at ksT/t = 0.01 for different angles tp with jq|a = 0.1. 
The result obtained from the Dirac-cone approximation is also 
shown (dashed line). Inset: Energy region around the van 
Hove singularity of the same curves. Right hand side: The 
imaginary part of the polarizability ImP' 1 ' (q, lj) as function 
of the energy u) at ksTjt — 0.01 for various wave vectors q 
defined in the text. 



only the subsequent integration over y has to be per- 
formed numerically. 

On the left hand side of Fig. 2, the imaginary part 
of the polarizability ImP^ (|q|, tp, ui) as function of the 
energy w is shown for different directions of the incom- 
ing wave vector q with |q|a = 0.1, where the usual 
parametrization in terms of the polar angle tp with q x — 
|q| cos tp is used. There is no apparent angle dependence 
except for the region around the van Hove singularity 
which is highlighted in the inset. The result obtained 
from the Dirac cone approximation is also shown (dashed 
line), which is given by 12 



ImR 



(1) 



0, Dirac 



i 



(14) 



For low energies, there is good agreement with the above 
formula, but especially for energies close to the van Hove 
singularity, u = 2t, strong deviations are seen which shall 
be discussed in the following in more detail. 



Expansion around the van Hove singularity 



where we defined 



hj(x,y;q x ,q v ) = \<pj(x,y;0,0)\ + \(j>j{x,y;q x ,q y )\ - \u>\ 



with 



(9) 



\<f> j (x,y;q X) q y )\ = 3 + 2(1 - y 2 ) cos(2^) (10) 



- 4^1 - y 2 ysm(2q y ) + j4(V 1 - x 2 cos(q x ) 



xsm(q x ))(y/l- y 2 cos(q y ) - ysm(q y )) 



1/2 



and 



F-(xva a)~- 1 F {x,y-q xi q v ) 



(11) 



with 



F j {x,y;q x ,q y ) = cos(2g x /3) + j2 y/l - y 2 (12) 
x (y/l - x 2 cos(2g x /3) - x8m(2q x /3)\ 

+ 2(2^1 - y 2 cosfe/3) + j [Vi - x 2 cos(fe/3) 

~xsin(q x /3) ) (y/l-y 2 cos{q y ) -ysm(q y )j . 

Above, we also introduced q x — 3q x a/2 and q y = 
y3q y a/2. Furthermore, the sum over x% is over all ze- 
ros which satisfy 



hj(xi,y;q x ,q y ) = 



(13) 



which can be written as a polynomial of fourth order. 
The zeros X{ can thus be obtained analytically such that 



The new feature compared to the Dirac cone approxi- 
mation comes from the region around the van Hove sin- 
gularity, located at the M-points of the Brillouin zone. 
For the Brillouin zone defined above, the M-points are lo- 
cated at M = 2tt/(3o)(1,0) and M± = n/(3a)(l,±y/3). 

In the following, we introduce the substitutions p x = 
3p x a/2 and p y = \/3p y a/2 and shall assume p x ,p y <C 
1 for p = k,q. Expanding around the M -point, the 
dispersion then simplifies to 



k 2 



"1 % 2 kff: 

and the band-overlap yields 
/ Mo (k,q); 
For the M±-point, we obtain 



4^/9 



1^1 



1 I 2 kj^ ky 



2ki 



and the band-overlap yields 



/_ f± (k,q) 



(q x /3±q y 



% \2 



(15) 
(16) 

(17) 

(18) 
(19) 

(20) 



For q y = (tp = 0), an analytical approximation sim- 
ilar to that presented in Ref. 29 is possible for the Mo- 
point expansion since the polynomial in the delta func- 
tion is quadratic. This yields a logarithmic divergence at 
wjw/t = ujm = 2 + <7 2 /2 which can be approximated by 
the following expression: 



lmPW> M °(q x ,i 



2sgn(w) Try/3 ( 2 
(2tt) 2 T \3a 
8A 2 



In 



9x 



In 



2 ~2 

18 

2A 2 g 
(Q - ui M ) 2 



(21) 
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where A denotes a suitable cutoff. Details on the calcu- 
lation are given in the appendix. For general q or for the 
M±-point expansion, we expect a similar behavior. 

In the inset of Fig. 2, the region around the van Hove 
singularity which is highlighted. For a general angle tp, all 
three M-points contribute and there is a prominent dou- 
ble or even triple peak structure. But for tp = tt/6, it/2, .. 
the overlap function of one Ai-point vanishes and the 
peaks merge. 



2. Behavior at large q 



lqla=0.O? , n=0.05f 



lqla-0.0? , n=0.5f 




0.014 0.015 
<B/f 



Let us now discuss the behavior for general q at /i = 0. 
For that we expand the energy dispersion |^k+q| around 
the points of high symmetry k = S = T, K, M and de- 
termine the q- vector for which \<j)^\ =0. To discuss the 
dielectric function at large wave- vectors |q| ~ 1/a, also 
local field effects have to be taken into account, 30 as was 
done in Ref. 31 . 

Expanding the dispersion around T = (0, 0) and K = 
27r/(3a)(l,l/ v / 3), we find that for q = K, \^< K \ = 0. 
The spectrum of ImP^ 1 ) thus starts at u> = 0. Expanding 
the dispersion around the Afo-point, the wave vectors 
qf° = 27r/(3a)(0,l/ v / 3) and qf° = 27r/(3o)(l, 2/\/3) 
yield |0* /o | = and the spectrum of ImP^ thus starts at 



t. For the M±-point expansion, we obtain 







for q = M T . 

On the right hand side of Fig. 2, the imaginary part of 
the polarizability is shown for the above wave vectors q. 
The curves for the Af±-point expansion yield the same 
curves as the ones for the A/o-point expansion and are 
not listed. For comparison, we also show the behavior 
for one of the vectors which define the Brillouin zone, 
q = t>i, which (for T = 0) is identical to the density of 
states by rescaling w w/2. 



B. Intraband transitions 

For finite chemical potential /i > 0, there are also in- 
traband transitions. The extension of Eq. (8) to finite 
chemical potential n and finite T is straightforward. The 
main difference is that the function of Eq. (9) now reads 

hj(x,y;q x ,q v ) = \<f>j(x,y;0,0)\ - \</>j(x,y;q x ,q y )\ ± |w| 



(22) 



This expression suggests that there might be differences 
to the Dirac cone approximation also for small uj. In fact, 
even for wave vectors |q| and chemical potential \x for 
which the Dirac cone approximation holds (e.g. |q|a 
0.1 . fi/t — 0.05), we find deviations from the Dirac cone 
result. 

Let us first start the discussion by summarizing the 
results coming from the Dirac cone approximation from 



FIG. 3: (color online): The imaginary part of the polarizabil- 
ity lmP^ (\q\,tp,uj) with |q|a = 0.01 as function of the energy 
uj for various angles tp and two chemical potentials y,/t = 0.05 
(left) and n/t = 0.5 (right) at k B T/t = 0.01. Also shown is 
the analytical result coming from the Dirac cone approxima- 
tion at zero temperature (dashed line). 



Ref. 16 for which we introduce the functions 



1 



/(|q|, w ) , — — 

G > (x)^xy / x 2 - 1 



(3t|q|a/2)*| 

- cosh _1 (a;) . 

- cos -1 (2) , 



X > 1 

\x\ < 1 



(23) 



Due to the finite chemical potential, the polarizability is 



now given by ImFjJ^ = ImP (1 ^ rac 



ImAP 



(i) 



with us® 



fi. Dirac 



and 



3i|q|a/2 the additional term reads 



ImAP, 



(i) 



fj,, Dirac 



(|q|,w)=/(|q|,a;)x 



G>(^)-G>(2^) , W <^A.<2 /1 . 



G>(^) 
G<(i 




-2p 



, UJ > A UJ < 2/1 - 

, uj < ujj/ Auj > uj® — 2/i 



q 



q 



,uj > uj^ Auj < u>~ + 
, otherwise 



2/.t 



(24) 



A distinct signature of non-interacting 2D electrons in 
graphene is a divergent behavior of the polarizability or 
charge susceptibility at the threshold for the excitation of 
electron- hole pairs at u>® , see Eq. (14). This divergence 
is also present for gated or doped graphene with /x > 
and has been usually attributed to the absence of curva- 
ture in the spectrum. But even in the regime where the 
Dirac cone approximation does not hold, i.e., curvature 
in form of trigonal warping has to be taken into account, 
we find a divergent behavior at oj d for q = q x . For arbi- 
trary direction, we find a peak splitting even in the Dirac 
cone regime. 

Let us discuss the polarizability using the full Brillouin 
zone in more detail. As stated above, PW(q, uj) is in- 
variant under rotation of it /Z - independent of the chem- 
ical potential and we also find the symmetry of Eq. (7). 
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FIG. 4: (color online): The imaginary part of the polarizabil- 
ity ImP^dql, tp, oS) with |q|a = 0.01 as function of the energy 
uj for various angles tp and two chemical potentials (x/t = 1 
(left) and fi/t = 1.5 (right) at k B T/t = 0.01. Also shown is 
the analytical result coming from the Dirac cone approxima- 
tion at zero temperature (dashed line). 



The numerical results moreover suggest that for moder- 
ate chemical potential /i < t the angle-dependent polar- 
izability can be described by a single function where the 
angle tp only enters as parameter. 

In Fig. 3, lmP^\q,tp,u) with |q|a = 0.01 as func- 
tion of the energy uj is shown for various angles ip and 
two chemical potentials n/t = 0.05 (left) and fi/t = 0.5 
(right) at ksT/t = 0.01. Only in the direction of q x , 
i.e., (p = 0, there is agreement with the analytical result 
of Eq. (24) coming from the Dirac cone approximation 
at zero temperature. Interestingly, this is also the case 
for a large chemical potential fi/t — 0.5, where trigonal 
warping effects should come into play. For arbitrary di- 
rection, a double-peak structure appears even for small 
chemical potential [ijt = 0.05 for which the Dirac cone 
approximation should hold. 

In Fig. 4, the same quantities are shown for larger 
chemical potentials [i/t = 1 (left) and \ijt = 1.5 (right). 
The differences to the analytical result coming from the 
Dirac cone approximation (dashed line) now become ap- 
parent also for q = q x (ip — 0). For /i/t — 1 they manifest 
themselves at lower energies uj < uj d and for [ijt = 1.5, 
a double-peak structure emerges. But interestingly, the 
divergence still occurs at uj m uj® in both cases. 

The above curves were obtained for ksT/t = 0.01, 
thus slightly larger than room temperature, but we have 
also investigated the effect of different T. We find that 
the curves for tp — are basically unaffected by tem- 
perature, but that for arbitrary direction the algebraic 
divergences seen for \ijt = 0.5 or /i/t = 1 are smeared 
out at larger temperature as it is the case for n/t = 0.05. 
On the contrary, the curves for \ijt = 0.05 and tp ^ 
develop the algebraic divergence for decreasing temper- 
ature. Generally, we can say that the algebraic diver- 
gences become broadened when the energy set by the 
temperature is much larger that maximal peak-splitting 



at tp = 7r/6. Nevertheless, the peak splitting in directions 
of lower symmetry prevails also at elevated temperatures. 



IV. REAL PART OF THE POLARIZ ABILITY 

The real part of the polarizability shall be obtained 
numerically via the Kramers-Kronig relation 



1 r 6t 

ReP (1) (q,w) = - / rf W 'lmP (1) (q,w')- 

7T Jo W' 



2uj' 



(25) 



The left hand side of Fig. 5 shows ReP*- 1 ' (|q| , tp, uj) for 
energies close to the van Hove singularity with |q|a = 0.1 
and = for different angles tp. As expected, there are 
strong deviations with respect to the result coming from 
the Dirac cone approximation and the functions become 
negative. This opens up the possibility of the emergence 
of an additional plasmon mode since the plasmon disper- 
sion in the RPA-approximation is given by the relation 



+ t; q P (1) (q,w) =0 



(26) 



where eoo denotes the effective dielectric constant in- 
cluding high-energy screening processes. Since the ex- 
periments in Ref. 26 were done on suspended graphene, 



we set e c 



1. For the Coulomb interaction, we set 



|q|u q = e 2 /2e = 90eVA« Whv F w 24ta. For |q|a = 0.1, 
ta 2 /vq — 0.004 never crosses one of the several curves 
which all tend to zero for larger energies. With the bare 
hopping amplitude t w 2.7eV, we do thus not find an 
additional pole in the RPA-susceptibility. 

Still, there is a renormalization of the hopping ampli- 
tude which comes from the wave function renormaliza- 
tion of the 7T electrons. Near the van Hove singularity, 
this renormalization will be large and the matrix element 
will be reduced. With a renormalization of t — > 2t/3, we 
would find an additional pole in the plasmon dispersion, 
consistent with experiments. 

We also expect deviations from the simple one-particle 
spectrum around uj — 2t, because expanding the effective 
screened Coulomb potential within the RPA approxima- 
tion, we have 



Im 



'1 



= Im 



e + ?vP (1) (q,w) 

-K/eo) 2 ImP«(q,w) . 



(27) 
(28) 



As can be seen from Eq. (21), there is a logarithmic 
divergence at uj]\i/t = 2 + (3q x a) 2 /8 even for q x — » 
since the prefactor from the band-overlap is canceled by 
v 2 

The right hand side of Fig. 5 shows the real part 
of the polarizability ReP^dql, tp,ui) with |q|a = 0.01 
and fi/t = 0.05 for different angles tp. For these pa- 
rameters, the Dirac cone approximation is supposed to 
hold, but strong deviations are seen for tp > as in the 
case of the imaginary part. This opens up the possi- 
bility of a modified plasmon dispersion as discussed in 
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\q\a=0. 1, \i=Q \q\a=0.01, p.=0.05t 




1.98 2 2.02 2.04 0.013 0.014 0.015 0.016 0.017 



m/t to/f 

FIG. 5: (color online): The real part of the polarizability 
ReP^ 1 - 1 (|q|, ip, oj) as function of the energy 10 for various an- 
gles cp at fcsT/t = 0.01. Left: For energies close to the van 
Hove singularity with [q|a = 0.1 and p — 0. Right: For low 
energies with |q|a = 0.01 and p/t = 0.05. Also shown is the 
analytical result coming from the Dirac cone approximation 
at zero temperature (dashed line). 



0) = ^(J>/t. Since the static value of the polarizabil- 
ity enters in the expression of the screened Coulomb po- 
tential, it is independent of the polar angle ip, consis- 
tent with group theory. Second, there are additional ze- 
ros of the real part of the RPA-diclectric function, i.e., 
to 2 e 00 /u q = 0.0004eoo (set e.g. f M » 2) crosses the 
curves at various energies different from the one asso- 
ciated with the Dirac cone approximation. Nevertheless, 
they lie in the region where ImP^ 1 ) is finite and the new 
plasmon modes are thus damped. 

For the present parameters, the undamped solution oc- 
curs at slightly larger energy compared to the Dirac cone 
approximation but is independent of the direction ip. For 
larger wave numbers |q|a w p/t, stronger deviations are 
seen, i.e., the plasmon dispersion depends on ip, is more 
strongly damped and eventually vanishes. 



V. CONCLUSIONS 



\q\a=0.01, \i=1.5t 
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FIG. 6: (color online): The real part of the polarizability 
ReP^ (|q|, ip, (jS) for |q|a = 0.01 as function of the energy lu 
for various angles ip at fcsT/i = 0.01. Left: For the chemical 
potential p/t = 1. Right: For the chemical potential p/t = 
1.5. Also shown is the analytical result coming from the Dirac 
cone approximation at zero temperature (dashed line). 



Ref. 28 . But for the present parameters, we do not find 
an additional zero in the RPA-dielectric function, i.e., 
to 2 e 00 /u q = O.OOfMeoc crosses all curves at the same en- 
ergy (choosing, e.g., the high-frequency dielectric con- 
stant of silicon = 2). Nevertheless, for larger wave 
numbers |q|a p/t, deviations are seen, i.e., the plas- 
mon dispersion is more strongly damped and eventually 
vanishes since the square- root singularity is smeared out. 

In Fig. 6, the real part of the polarizability is shown for 
two large chemical potentials p/t — 1 (left) and p/t = 1.5 
(right). As for the imaginary part, large deviations com- 
pared to the results coming from the Dirac cone approx- 
imation (dashed line) are seen. First, the static value 



pW(|q|, 



ip,L) 



0) is larger than ta 2 R.eP^ irac (|q|, u> 



We discussed the polarizability of graphene using the 
full band structure of the 7r-electrons. We especially 
focused on the features around the van Hove singular- 
ity since they might be responsible for the newly found 
plasmon dispersion. We find that there are no plasmon 
modes coming from the van Hove singularity within the 
RPA-approximation with the bare hopping amplitude t. 
But with a renormalization of t — > 2t/3, there are ad- 
ditional plasmon modes, consistent with experiment. 26 
We also find a logarithmic divergence of the imaginary 
part of the effective Coulomb interaction. This will lead 
to prominent electron-hole interactions as was recently 
found in Ref. 27 . 

We also looked at the intraband contribution to the po- 
larizability. For q in the T — M-direction, we find basic 
agreement with the results of the Dirac cone approxima- 
tion even for rather large chemical potential p s» t/2, 
i.e., when corrections to the linear Dirac spectrum are 
large. For arbitrary direction of the incoming wave vec- 
tor q, we surprisingly found strong deviations from the 
results coming from the Dirac cone approximation where 
a double-peak structure emerges. The peak splitting oc- 
curs for all values of p > and is a direct consequence 
of the energy dispersion. For fixed q and p, it is largest 
for ip = 7r/2 and tends to zero for |q|, p, — > 0. As a conse- 
quence, the plasmon dispersion is more strongly damped 
for |q|a ~ p/t and eventually vanishes at larger tempera- 
tures since the square-root singularities are smeared out. 
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VII. APPENDIX 

For q y = and using the saddle-point approximation of 
Eq. (16), an analytical calculation of the imaginary part 
of the polarizability around the M -point is possible. The 
integral is preformed in polar coordinates and with x = 
cos ip, we have to solve the following quadratic equation 
with respect to x: 

\<fl?°\ + 0k+,. | - w = 4fcV ± 2kq x x + e = (29) 

with e = 2(1 — k 2 — Cb/2 + q 2 ), Co = \ui/t\ and q x = 3q x a/2. 
The integral over <p eliminates the delta function and 
yields 

-^'-fi^l't (30) 



Z+(<L,u>) + X_(<? x ,d>) 



(31) 



where we included the spin-degeneracy g s — 2 and de- 
fined the following integrals: 



l±(q x ,ui) = / 
Jx 



dkI±(k;q x ,oo) , (32) 



with 

I±(k;q x ,Co) = 



(33) 



and the integration domains 

' [fc min ,A] ; u <2 + 5q 2 x /8 



[fc_,A] ;w> 2 + 5^/8 



(34) 



and 



V-(q x ,u}) = < 



f [fcnin.A] ; u <2 + q 2 /2 

[fc mi „,fc_]U [fc+,A] ; 2 + ql/2<ij 
<2 + hq 2 x /8 
{ [k+,A] ;Cb>2 + 5q 2 x /8 



(35) 

We further defined £ 2 = -<g/4-l+w/2, k^ in = q 2 /8-£ 2 , 
k± = \£±q x /2\, and A denotes a suitable cutoff. 
The indefinite integral has the following solution: 



/ 



dkI±(k;q x ,oj) = - In 



±q x 



(36) 



+ \/8(fc 2 " k 2 min ) + \ Hk 2 e) ± 2q x J8(k 2 klJ 



There is a logarithmic singularity at % = 2 + q 2 /2. 
Expanding the above result around % yields the simple 
expression 
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